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Abstract 

Zeilberger |12] proved the Refined Alternating Sign Matrix Theorem, which gives 
a product formula, first conjectured by Mills, Robbins and Rumsey [9^, for the num- 
ber of alternating sign matrices with given top row. Stroganov ^Oj proved an explicit 
formula for the number of alternating sign matrices with given top and bottom rows. 
Fischer and Romik [7] considered a different kind of "doubly-refined enumeration" 
where one counts alternating sign matrices with given top two rows, and obtained 
partial results on this enumeration. In this paper we continue the study of the 
doubly-refined enumeration with respect to the top two rows, and use Stroganov's 
formula to prove an explicit formula for these doubly-refined enumeration numbers. 

1 Introduction 

An alternating sign matrix (ASM) of order n is an n x n matrix with entries in 
{0, —1, 1} such that in every row and every column, the sum of the entries is 1 and 
the non-zero numbers appear with alternating signs. See Fig. for an example. 
A monotone triangle of order n is a triangular array {ti^j)i<i<n,i<j<i of integers 
satisfying the inequalities 

ti,j+i {2<i<n,l<j <i-l). 

A monotone triangle of order n is called complete if its bottom row has the numbers 
(l,2,...,n). See Fig. Qb). 
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Figure 1: An ASM of order 5 and the corresponding complete monotone triangle. 



For integers ki < k2 < ■ ■ ■ < k„ 



monotone triangles with bottom row (/l]^ , . . . , /c- 

An = a„(l,2, . . . 
An,k = a„_i(l,2,...,A:,...,n), 
An,i,j = a„_2(l,...,i,...,j,...,n), 



denote by an{ki, k2, ■ ■ ■ , k„ 
Define 



the number of 



(1 < A; < n) 

(1 < i < j < n). 



where the notation s in a list of numbers indicates that s is omitted from the list. 
For notational convenience later on, take A^^k as if A; ^ [1, n]. 

Since ASMs of order n are well-known to be in bijection with complete monotone 
triangles of order n, the number A^ is the total number of ASMs of order n. From 
obvious properties of the bijection, it follows that An^k is the number of ASMs 
of order n such that the unique 1 in the first row is in position k. The numbers 
n^k s-re called the refined enumeration numbers for ASMs. Zeilberger jll] 112] 
proved the celebrated alternating sign matrix theorem and refined alternating sign 
matrix theorem, which state, respectively, that 
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In [7j the numbers (^n,ij)n,i,j were studied, in an attempt to extend the work of 
Zeilberger to a doubly-refined enumeration of alternating sign matrices based 
on their first two rows. We refer to these numbers as the doubly- refined enumeration 
numbers. The connection to enumeration of alternating sign matrices (discussed at 
length in [7]) is illustrated in Fig.[2| which shows the possible configurations for the 
first two rows of an ASM. For fixed n, these configurations are indexed by a triple 
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Figure 2: The possible configurations for tlie first two rows of an ASM: (a) i < k < j. (b) 
i = k < j. (c) i < k = j. In all three cases the number of ASMs with the given first two 
rows is Ar 



k) such that l<i<k<j<n and i < j, and for each such triple the number 
of ASMs having the given first two rows is independent of k and in fact is equal to 
A ■ ■ 

The results of [7] gave only partial information on the A„^jj 's, namely a system of 
linear equations satisfied by {An^ij)ij for each n. These equations, along with some 
other information known from simple considerations, were conjectured to determine 
the An^ij''s uniquely and so to allow to express the An^ij^s as ratios of determinants 
and to compute them efficiently on a computer for reasonably large values of n. A 
rather complicated explicit formula not involving determinants was also conjectured 
based on numerical evidence. 

Our main result is a new and simpler explicit formula for A^^ij. We prove: 

Theorem 1. Let 

For each l<i<j<nwe have 

n-j p 



"■-J P / \ 

p=0 q=0 



(4) 



Our proof of Theorem [T] uses the well-known connection between alternating sign 
matrices and the square ice model, and builds on previous results and techniques 
of Izergin-Korepin [3 , Kuperberg I^Sj and Stroganov [lOj. Following our discovery 
of Q, an alternative derivation of the same formula using the monotone triangle 
techniques developed in |H [5l [7] was recently found by Fischer [6] . 

The following conjecture holds empirically but does not follow from our methods. 
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I 1 t 1 t 1 

(type 1) (type 2) (type 3) (type 4) (type 5) (type 6) 
w{v) = a w{v) = a w{v) = b w{v) = b w{v) = c w{v) = c 
m = 171 = 171 = m = m = —1 m = 1 

Figure 3: The six types of vertices in a square ice configuration, the weight w attached 
to each of them and the entry m assigned in the corresponding ASM under the bijection 
translating square ice configurations to ASM's. 



Conjecture 2. For each n > 1 let (An^ij) be the extended doubly-refined 

enumeration numbers defined in '^jj. Then d4| holds also for the extended numbers. 
In other words, we have 

^n,i,j = Y.Y.^-iyr]Xnii + q,j+p), (1 <i,j <n). 

2 Square ice and the partition function 

We consider square ice (or six- vertex model) configurations on an n x n square 
lattice 

L„ = {1, . . . ,n} X {1, . . . ,n} 

satisfying the so-called domain wall boundary conditions. In a square ice con- 
figuration the edges of the lattice are oriented so that each vertex in the lattice 
has two incoming edges and two outgoing edges, giving six possibilities, shown in 
Fig. |3j The boundary conditions are that edges entering the lattice from the left 
and right sides of the square point inwards, whereas edges adjacent to the top and 
bottom sides point outwards. Such configurations are in bijection with alternating 
sign matrices of order n. The bijection maps each of the six types of square ice 
vertices to either a-|-l,a— loraO; see Figs. [3]and|4| 

Fix a real value r] called the crossing parameter. Given a square ice configu- 
ration C, and given a vertex v in the lattice with an associated (real-valued) row 
parameter x and column parameter y, according to the type of the vertex v in 
the configuration C we associate with it a weight w{v) = wc{v) equal to either of 
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Figure 4: The square ice configuration corresponding to the ASM in Figure ITVa) . 



the three quantities a, b or c, defined by 

sin(ry/2 + x — y) 
sin T] ' 
sm{rj/2 — X + y) 
sin r] ' 

1. 

Vertices of type 1 and 2 are assigned the weight a, those of type 3 and 4 are 
assigned the weight 6, and types 5 and 6 get the weight c. The weight of a square 
ice configuration C is then defined as the product of the vertex weights over all the 
vertices in the lattice, namely 

w{C) = Y[ wciv). 

Denote the set of square ice configurations by C„. The square ice partition func- 
tion is defined as the sum of the weights over all configurations and is denoted by 
Zn, i.e., 

cec„ cec 

Here, we associate with a vertex in row i and column j a row parameter Xi and a 
column parameter yj. Thus the partition function Zn is a function of the crossing 
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parameter rj and of the variables xi, . . . , Xn,yi, ■ ■ ■ ,yn, called the spectral pa- 
rameters. Regarding 77 as fixed, we occasionally write Zn{xi, . . . , Xn', yi, ■ ■ ■ , Vn) 
to emphasize the dependence on the spectral parameters or to substitute specific 
values. 

An important fact that we will use is that Zn is a symmetric function in the row 
parameters xi, . . . ,Xn and is also symmetric in the column parameters yi, . . . ,yn- 
For the remainder of the paper we take r] = 27r/3. For this particular value of the 
crossing parameter, Stroganov |10! showed that Zn is actually symmetric in the 
union of the Xj's and y/s, but we will not need this fact. 

For our proof of Theorem [T] we will consider two specializations of the partition 
function, namely 



where u and v are parameters. By the symmetry of Zn in the row parameters, 
we get that Si = 82- This will enable us to relate the doubly-refined enumeration 
numbers to a different kind of doubly-refined enumeration studied by Stroganov, 
where one enumerates ASM based on the first and last row instead of the first two 
rows. 



For n >1 and 1 < i, j < n, let Bn.ij denote the number of ASM's of order n with a 
1 in positions (1, i) and (n, j). These numbers enumerate ASM's based on their first 
and last row and so are also sometimes referred to as the doubly-refined enumera- 
tion numbers. They were originally considered by Mills, Robbins and Rumsey [9]. 
Stroganov proved the following formula expressing -Bn,i,j in terms of the (singly-) 
refined enumeration numbers An^k (see also [2] where some generalizations of this 
result are proved). 

Theorem 3 (Stroganov [TO]). We have 

Bn,i+l.j+l~Bn,iJ =Yn{'i':j) •= ^, j' (^n,j+l ~^n,j) ~l"^»i-l,i(^n,j+l ~^r»,j ) ■ 

A^-l y / 

Note that Yn{i,j) = Xn{i,n — j) where Xn{s,t) is defined in Q. By summation 
of the differences it follows also that 

min(i,j) — 1 



Si 

S2 



Zn{u, 0, . . . , 0, 0, . . . , 0) and 
Zn{u,v,0,...,0;0,...,0) 



3 Evaluation of S 
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Stroganov's proof involved evaluating the expression denoted above by Si in two 
ways. For our purposes, we only need one of them, whose derivation we include for 
completeness. Let 

^(.) = -^sm(./3 + x), t = tiu) = ^^, . = = 

Now note that when summing over all square ice configurations to evaluate Si = 
Zn{u, 0, . . . , 0, f ; 0, . . . , 0), the weight of each configuration depends only on the po- 
sitions of the 1 's in the first and last rows of the ASM corresponding to the configu- 
ration, since all vertices in rows 2 through n—1 contribute a factor of 1. If the ASM 
has I's in positions i and j respectively in the first and last row, using the transla- 
tion between ASM's and square ice configurations the weight of the corresponding 
configuration is easily seen to be 

This implies that 

n 

— 1 " 

= {cp{vM-u)y~ Y.Bn,ijf-'s^-'. (5) 

4 Evaluation of ^'2 

To evaluate ^2, note again that for this substitution of spectral parameters the 
weight of a configuration is only dependent on the state of the first two rows of the 
square ice configuration, or equivalently of the corresponding ASM. As explained 
above, this state can be indexed by parameters k satisfying l<i<k<j<n 
and i < j, which correspond to the three numbers in the top two rows of the complete 
monotone triangle that corresponds to the ASM. The weight of the corresponding 
square ice configuration having given parameters i,j,k is then easily computed to 
be 

'^{u)''-'^ip{-u)''-''^{vy+^-^-^ip{-v)''-^+^-'-^ i<k<j, 

w{c) = I ip{u)''-^ip{-uy'-''(p{vy-^(p{-v)''-^+^ i = k<j, 

^Lp{u)^^^Lp{-uy^^^Lp{vyLp{-v)'^~'^~'^ i < k = j. 

Now sum this over all configurations, and use the "120-degree triangle" identity 

(p{xf + (p{-xf - (p{x)(p{y) = 1, 
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to get that 
S2 = 



l<i<j<n \ 



+99(n)^-V(-?i)""V(^')V(-^)""'"' 

k=i+l / 



l<i<j<n \ 
l<i<j<n \ 

i-1 



+ (1 + 5^-5) Y t^-^ S'''^+^-'-^ . 

k=i+l / 



The inner summation on /c is a finite geometric series that sums to 

~ ts-1 



J ^ n-i-1 j-i„ri-j 

E^k-l^n-j+k-i-1 _ ^ ^ ^ 

t<i - 1 

k=i+l 

so that, after some further simple algebraic simplifications, we obtain 

{tS - 1)52 _ A . .(fi^n-j+l _.i n-j n~j+l 

+ t^s""-' + - t^-h""-' ) . (6) 

5 Completion of the proof 

Now knowing that Si = S2 we can write 

{ts - l)Si {ts - 1)82 



We can equate the coefficients of fs^ on both sides (since these functions are linearly 
independent, a fact that is equivalent to [T, Ex. 7.1.9, p. 231]). In conjunction with 
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^ and ^ this translates to the identity 

'^nihj) = Bn,i+l,j+l — Bn,iJ 
— ^n,j+l,n+l— J ~l~ ^ra,i,n— j -^n,i,n+l—j -^n,n—j,i -^n,n+l—j,i+l ~\~ -^n,n—j,i+li (7) 

which holds for all i, j if we adopt the convention that A^^p^q and Bn^p,q are outside 
the respective ranges of their definitions. 

It remains to solve this system of equations in the An^i^s. First we reformulate 
it slightly for convenience by replacing j hy n — j. This gives 

^nihj) = Ynih n — j) = Bn,i+l,n-j+l — Bn,i,n-j 

This simplifies even further when one observes that the last 3 terms on the right- 
hand side are in the range of parameters 1 < i < j < n which interests us. So we 
have 

Xn{i,j) = + l + An,i,j - An,i,j+l, {'i- < i < j < n) . 

Rewriting this in the form 

it can now be easily solved by iteration (or more formally by reverse induction on 
j), as follows: 

= Xn{i,j) + Xn{hj + 1) - Xn{i + 1, j + 1) 

+ (^n,j,j+2 — ^n,i+lj+2) " (^n,j+l,j+2 " ^n,i+2,jr'+2) 
= Xn{i,j) + Xn{iJ + 1) - Xnii + 1, j + 1) + " 2An,i+lJ+2 + An,i+2,j+2 

= Xn{i,j) + {Xn{i,j + 1) - Xn{i + + 1)) 

+ {Xn{i,j + 2)-2Xn{i + l,j + 2)+Xn{i + 2,j + 2)) 

+^n,«j+3 — 'iAn,i+l,j+3 + 3^n,i+2,j+3 — ^n,j+3,j+3 
A:=0 ^ ^ p=0 q=0 ^'^^ 

In the last expression, the first summation over k corresponds exactly to the case 
p = n — j oi the second sum, so we can shorten this to 

= EE {\-^fXn{i + q^J +P) 

p=0 q=0 

which was the claim of Theorem [TJ □ 
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